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Abstract 

The stochastic quantization of the fermion field is performed starting from Dirac 
equations. The statistical properties of stochastic terms in Langevin equations are 
described by explicit formulae of a Markov process. The interaction of the field is 
introduced as correlation of the stochastic terms. In the long time limit free fermions 
disappear and proper combinations of field components propagate as a scalar boson 
field. The existence and uniqueness of the long time limit is proved in the first order 
approximation of stochastic Liouville equation. 
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1 Introduction 



The method of stochastic quantization [Q, 0] starts from Euchdean field equations amended 
by addition of a random stochastic variable. Thus the change of the field during an 
infinitesimal increment of stochastic time t which is considered to be an auxiliary variable 
is described by field equation and stochastic term. The massive scalar field is described 
e.g. by 

t) = i-d' + m2)0(x, t) + t) (1.1) 

where x are coordinates in Euclidean space. ( p..l| ) can be also regarded as Langevin 
equation with drift term equal to Klein-Gordon operator. The statistical properties of 



the fluctuation term ^{x,t) are given. One particular solution of ( |1 . 1| ) is of no special 
interest but the mean value of the correlation {(j){x,t),(j){x ,t)) averaged over all possible 
fluctuations is equal in the long time limit (i.e. t — > oo ) to the usual result of path 
integral quantization of the Euclidean field theory under consideration. One has of course 
to prove that the long time limit exists and it is unique. 

In the standard approach the fluctuation term corresponds to white noise - i.e. it is a 
random Gaussian variable with autocorrelation equal to 

{^{x,t)^{x',t')) = c6{x -x')6{t-t') (1.2) 

The finite width of the Gaussian distribution of ^ in t variable in (|1.1|) leads to Pauli-Villars 
regularization of propagators 0. 

Stochastic quantization of boson fields is the easier one. Stochastic quantization of 
fermion fields is more elaborate and there are two approaches. The method proposed 
by Breit, Gupta and Zaks 0] starts also from Klein-Gordon operator as in (|1.1|) and the 
fermion propagator results from special statistical properties of the fluctuation term. In 
this approach 

d 

—i){x, t) = (-9' + m^)^{x, t) + ri{x, t), 

d - - 

—ip (x, t) = V {x, t)(-<9^ + m'^) + f] (x, t) 
at 

and 

{ri{x,t)f] {x' ,t')) = c{-i^ +m)S{x-x')5{t-t'). (1.3) 

Thus quantization of boson and fermion fields is put on the same footing and the stochas- 
tic time has in both two cases the same dimension [t] = m~^. This is convenient for 
perturbative description of interacting fermion fields. 

The second method proposed by Fukai et al. p starts from Dirac operators 

d 

—i){x,t) = {i^ -m)ij{x,t) +i{x,t), 



;i-4) 



d - - 

—ip {x, t)=i) (x, t){-i^ - m) + ^ (x, t) 



and the stochastic time t has the real time dimension. It seems to be quite difficult to 
implement this method in the case of interacting fermion fields within the framework of 
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standard perturbation approach. However, it allows to quantize fermions in such a way 
that correlations of proper fermion field combinations are equal to autocorrelations of a 
boson field. Such procedure is the subject of this paper. 

The idea is to introduce the interaction between fermion fields via correlation of fluc- 
tuation terms. Direct consequence of simple correlations of the type ( |1.2|) is locality of 
resulting theory. As the correlations of the fluctuation terms are related to those of the 
field components, more elaborate probability distribution of stochastic terms can lead, in 
addition to ensuring the locality, also to dynamical consequences. To achieve this it is 
necessary to use one multivariate distribution governing behavior of all fluctuation terms 
in Langevin equations at given stochastic time t. The sequence of these distributions will 
form a Markov process. It is evident that random variables in ( |1.4| ) have to be of anticom- 
muting nature. Markov processes of Grassmann numbers are constructed and studied in 
the Section |^ of the paper. Lowest order propagator based on Langevin equations with 
correlated fluctuation terms is calculated in the Section |[ The existence and the unique- 
ness of the long time limit of the studied system is proved in flrst order approximation in 
the Section ^. The Section | of the paper is devoted to conclusions and discussion. 



2 Markov chains of Grassmann numbers 

Classical Markov processes Q are based on the notion of probability distribution. However, 
there is no such concept as probability measure in Grassmann number calculus. Hence 
one has to adopt a more general approach outlined e.g. in and also used in IQ. It is 
possible to use as a probability distribution of Grassmann variables ^, ^ (henceforth it 
will be dealt always with pairs of Grassmann numbers) a function which can be generally 
written as 

P{^, a = Co + + C2^ + cs^ (2.1) 

and such that the integral 



gives correct expectation values (/(^,0) function /. It is reasonable to assume 

that P commutes with any Grassmann number. The conventions of Grassmann variables 
integration and derivations (left-derivatives) are the same as those used in 0, 

d^ = 0, Udi = i, I d^ =0, =i (2.3) 



so that e.g. 

" iUidi =1. 

|(ee ) = e , |-(a ) = --^{i = ■ (2.4) 

Classical Markov process is given by the probability distribution P{yi,ti) of stochastic 
variable yi at certain time ti ( with the abbreviation yi = y{ti) ) and the transition 

-'^An excellent introduction into Markov processes and stochastic differential equations can be found 
in§. 
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probability T{y2, t2 \ yi,ti) which determines the probabihty distribution of y at any time 

t2>ti 

P{y2, t2) = I T{y2, t2 I 1/1, ti)P(z/i, ti) dyi (2.5) 
The transition probabihty satisfies Chapman-Kolmogorov equation for any > ^2 > 

7^(2/3, ^3 I 2/1, h) = J T{y3, ts I 2/2, t2)T{y2, ^2 | 2/i, h) dy2 (2.6) 
and it is normahzed 

jT{y2,t2\yuh)dy2 = l. (2.7) 

One of the most simple stochastic processes one can use is the Orenstein-Uhlenbeck pro- 
cess. The natural step therefore will be to construct Grassmann variable analog to this 
process which is the unique stationary Gaussian Markov process and it is defined by 

PiVi) = e-^, (2.8) 



T{y2M I yiM) = I , ^ == exp 



(2/2 - Vie 



-Mt\2 



2a Jl - e 



-2Mt\ 



(2.9) 



where t = t2 — ti. 

The multivariate distributions are needed but for simplicity in the following the func- 
tions of only one pair of Grassmann variables C, and will be used. All necessary 
mathematical properties can be demonstrated by means of such functions. One can take 
as Grassmann Gaussian probability distribution 

p(e,o = -i-ee (2.10) 

and readily check that the normalization condition holds 

/ =1. (2.11) 

If one tries as transition probability distribution 

^2,^2 I = -1 - ^2^2 - 6^^'^'^-'^%^, - ^"26] (2.12) 

then the following equations hold 

/ T(6, ^2 I ei, QPi^i, Q d^.di, = P(e2, ^2), (2.13) 

/ ^(6, ^3 I 6, ^2)^(6, ^2 I 6, Q di2dl2 = T(e3, ^3 I 6, (2.14) 

provided that ^3 > ^2 > , and finally 

/T(6,e2iei,ei)«2 = i- (2.15) 



They correspond to (|2.5|J2.6| ) and (|2.7| ) respectively. However, there is one significant 
difference - the exponent in ( p.l2| ) can be also positive in contrast to (|2.9| ) where the 
positive exponent in t dependence would lead to complex T values. 
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The analogy of ( ^.10| , pTT^ with the classical process can be extended even further. 
The equation 



d d - d'^T 



(2.16) 



is the forward Kolmogorov equation for Grassmann variable Markov process ( p.l2|) , while 
the forward Kolmogorov equation of the classical process (|2.8| , p.9|) reads 



dT 



= -M 



dy2 



The backward Kolmogorov equation for ( |2.12|) has no analogy in corresponding classical 
equation. The full analogy is restored for 



T (^2, ^2 I ^1; 



^±2A/r 



(^2 



6). 



:2.i7) 



This transition probability distribution is the full analog to that of classical process, 
including both the forward and backward Kolmogorov equation. In contrast with classical 
Orenstein-Uhlenbeck process, which is the unique stationary Gaussian Markov process, 
here one has more stationary processes which can be called Gaussian. This is only natural 
as one deals with truncated Taylor series. The Kolmogorov equations of the process of 
anticommuting variables ( p. IT] ) have the same asymptotics for r — > ±cxd as the classical 
process. In the following the analog of simpler ( p. 121 ) will be used. 

The stochastic variables which will enter the Langevin equations will be represented 
by quartets of functions which ascribe to each point in Euclidean space-time x and to 
each value of stochastic time t four random Grassmann numbers. To achieve the desired 
result - coupling of four fermions - it turns out that the probability distributions of these 
variables have to be functions of two pairs of such stochastic variables , ri^{x,t) 

and ^^{x',t), f]^{x',t), where /i = 0,1,2,3). These pairs will be defined on a product 
of two Euclidean spaces x x x' and their distribution in this space will be described by 



(5-function 5{x 



X 



Further one has to suppose that these "quartets of quartets" have 



nonzero product of all their components - i.e. all 16 Grassmann numbers are different. 
Thus the probability distribution in x' becomes for any x 

V{i{x\ t), e(x , t), r/(a;', t), n{x\ t)) 1= P{^{x, t), ^(x, t), r]{x, t),f]{x, t)) 6{x - x') (2.18) 

and the transition probability distribution as a function of y' will be 

^ t2),ay', h), Viy', h), f]{y', h) I i{x, h),^{x, h), r,{x, ti), n{x, h)) \y = (2.19) 
T{$,iy, t2),^{y, t2),r]{y, t2),f]{y, ts) I ^{x, ti), ^(x, ti), ri{x, ti), f]{x, ti))5{y - y') 



for any x, y. Both in (|2.18| ) and (|2.19|) the greek indices are omitted. In rigorous approach 
the probability distributions in Euclidean space in ( |2.18D and (|2.19|) should be e.g. of the 
type listed in (|2.8|). One can consider the 5-functions in ( p.l8| , p.l9D to result from the 
limit — *■ in (|2.8|) . Simultaneously the 0"^ parameter disappears. In fact one allows in 
( p. ID as c-coefficients not only c-functions but also c-distributions. 

As the random Grassmann variable distribution P{^{x,t),^{x,t),ri{x,t),f]{x,t)) in (|2.18|) 
at given stochastic time ti the following function will be taken as ansatz 



UvlUCUCUvl+ JliUvlUCUCUvl) (2.20) 



(7=0 



u=0 



i,j=Q M=o 



p=0 



<T = 



u=0 
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where the greek indices of ^, ^, and fj on the left hand side of ( p^.20[ ) are omitted and the 
X dependence is omitted on both sides. The stochastic time dependence is indicated by 
upper right index, e.g. f]{ti) = t]^. Transition probabihty distribution on right hand side 
of ( p.l9| ) written in the same notation reads 



where 



+ k E (-1)^-^^- n n c n n ■ + (2.21) 

jj=0 M=o p=o (T=0 y=0 

+ k E (-1)^-^^- fivi li c n e n v'. ■ viij , 

ij=0 ^l=0 p=0 o-=0 iy=0 



k = e^*^" , T = t2 - h , t2 > h. 



Distributions (p.l8| ) and ( |2.19| ) obey the relations ( |2.13| ) - ( |2.15| ) amended by correspond- 



ing space-time integrations preceeding the Grassmann variable integrations. This stochas- 
tic process leads to correlations 

{^p{Xi,ti) r]p{x'i,ti) f]^{x2, ^a{x'2, t2)) = Ce 6{xi - x' i) 6{X2 - x' 2) ^pa C^^'^^*^'^^^ 

(2.22) 

where ti < ^2- The constant = 1, [c^] = is introduced to adjust dimensions. 
As the stochastic variables will disappear after integrations in all applications in Sections 
3 and 4 it is simpler to consider them to be dimensionless and to adjust dimensions in 
Langevin equations by another constant r^,. This constant can be also used to readjust 
the numerical value of so that it remains equal to 1 even after change of units. 

All expectation values of single Grassmann variables or incomplete pairs rj oi fjE, are 
are equal to zero, e.g. 

{Ipi^iM)) = {Vp{xi,ti)) = {f]p{xi,ti)ri^{xi,ti)) = 

(2.23) 

= {^p{xi,ti)Vu{xi,ti)f]p{xi,ti)) = etc. = 0. 

Also if there is not complete pair rj^ or fjC, at any time t the correlation is equal to zero, 
e.g. 

{^p{xuh)r]^{x2,t2)f],{x3,ts)Ux3,t3)) = (2.24) 
for ti < t2 < ^3- The expectation values are zero also in case of the combinations like 

( ^pi^U h) r]p{Xi, ti) ^^{X2, t2) vAx2, t2) ) = 

(2.25) 

= {^piXi,ti)f]^{xi,ti)^^{x2,t2)f]^{x2,t2)) = 

assuming ti < t2. Also the mean value 

{Ip{.xi,ti)ip{x2,t2)f]^{x'i,h)'q„{xi,ti)) = (2.26) 
whenever the times ti, t2, and are not equal each other. 
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The higher moments are reduced to the product of the lower ones, e.g. 

{^piXuti)Vf^iXutl)Vuix2, t2)ia{.X2, t2)^^ix3, ts)ri^{x3, h)f]s{.Xi, ti)ip{Xi, U)) = 

{VpiXl,ti)^^{Xi,ti)^^{x2, t2)Vaix2, t2)^aix3, h)ri^{x'i, h)f]s{xi, t4)ip{xi, U)) (2.27) 
= {Vpi^u ti)^f,ixi,ti)^^{x2, t2)r]^ix2, t2)){^^ixs, t3)r]^{xs, h)f]s{x4, t4)^^{x4, U)) 

provided that ti < )f;2 < ^3 < ^4 • 



3 Coupled fermions 

The subject of the interest will be a fermion field i/j with no other interaction except the 
one represented by Markov process ( p.l8| ),( p.l9|) . The Langevin equations in this case are 



d 



(3.1) 

d - - - _ 

—ip^ = {i^ + m)^„ + (^„ + 7]^) 



where 

^ =7'^^. = 7^^ (3.2) 
and 7^ are Dirac matrices in Euclidean space, satisfying 

{7^7^} = -25,.. (3.3) 

The is a constant with dimension [m^/^]. The stochastic variables are dimensionless as 
discussed in Sec. 2. 

To solve (13.11) it is useful to introduce retarded Green functions 



G{x,t) = Q{t) J -^exp[-{^ +m)t + tpx] 



which satisfy 
and 



(2vr) 



j^-^ exp[-(|^ + m)t - ipx] 
G{x, t) = G {x,t) = ioi t < 0, 

d 

— Ga-y{x,t) - ~ m)af3G/3^{x,t) = 5a'r5'^{x) 5{t) 



d - 

— Ga-y{x,t) - G a/six, t){-i(^ - m) (3^ = 5a^5'^{x) 5{t) 



(3.4) 



(3.5) 



Solutions to ( p.l|) can be written as convolutions 







ijja{x,t) =ro / dr d x' Gap{x - x',t - t){^p{x',t) +'r]p{x',r)) 



(3.6) 

ilj^{x,t) =rD I dr I d^x' {L{x',t) + f] {x',T))Gpa{x - x',t - t) 







The nature of fluctuation terms in (|3.1| ) implies that the simplest nonzero correlation 
function one can calculate (see( p.22| )) is 



{'lpa{Xl,ti)llj°'{x2,t2)lp/3{yi,t3)tlj'^{y2,U)) = r'^^J dTiJdT2jdT3jdT4j d^ij d^J y y u, y 2 
{{^„{Xi,Ti)+f]^{x'i,Ti)) G^a (Xi - Xi , 1 1 - Ti ) G^p (^2 - X2 , )f:2 - ^"2 ) (^p , ) + ?7p (x'a , Ta ) ) 

i^Mi > ^3) + ^p , ^3) ) Gpf3 {yi - y[ , ts - rg) {y2 -y'2,U-T4) (^^ {y'2,r4)+r]^{y'2,T4))) 

= rl[dTi[dT2!dn!dTJd%td%tdyJd%'2 (3.7) 
Jo Jo Jo Jo J J J J 

Gaa (si - x'l , ti - Ti ) Gap {x2 - x'2 ,^2-^-2) G pp {yi - y\ , - Tz)Gfiy (?/2 - y 2 , ^4 - T4) 
( (^^ (x'l , Ti ) + r/^ (x'l , Ti ) ) (^p (X2 , r2 ) + ?7p (x^ , r2 ) ) X 

X {liivi^ ^3)+^^(i/i, r3)){iv{y'2^ u)+Vu{y2, u))) 

where (...) represents integration over Grassmann variables. Due to the properties of 
stochastic variables ( p.25| , p.26D , integral ( p^) over region ti x ^2 x ^3 x ^4 is equal to 
zero. To stay on the manifold where the mean value (...) 7^ 0, it is necessary to impose 
condition that the time integrations are carried out only for ti = t2 < = t4^. Hence 
( ^?7D will be first integrated with the weight 

CrS{Ti-T2)S{T3-Ti) , (3.8) 

where = 1, [c^] = m~^. The numerical value of Cr can be again adjusted to 1 by means 
of r^, constant as in all expressions leading to measurable quantities one has always the 
product Cr Ce which is dimensionless. The condition ( |3.8| ) means, that the product 
ipai^'^ will be always taken in the same stochastic time instant. From ( |3.7|) one arrives at 



{^Ij^r^pi;^) = 2crCEri J^dnJ^dTsJd^x'.J d'y\J 



d pi f d p2 f d f d 



{2n)y (27r)V (27r)V {2n)^ 



6 (ti - Ti) exp[-(jf^^ + m)(ti - ri) - ipix^ 

e (t2 - n) exp[-(|^2 + "^)(^2 - T\) + ^P2X2]ap exp[i(pi - P2)x'l] 

(^3 - T-3) exp[-(;^3 + m)(t3 - rg) - ip-iy^pp 

e (^4 - Tg) exp[-(;^4 + m)(t4 - Tg) + ip^y-^du exp[z(p3 - PA)y ^ 

(5,,, exp[±M I rg-Ti I] (3.9) 

carrying out also the integrations over 5-functions 

5(x'i -x'2)^{y\ -y'2) 
which come from ( p.22| ). Further integrations lead to 



i^ari^pi^') = ^rfdnfdrJ ^/ ^ e(Ti - r0e(T3 - ra) 
Jo Jo J [271 ) J {2tt) 

exp[-2(|^^ + m)(Ti - ri)]^pexp[-2(jf^3 + m)(T3 - rg)]^^ 

exp[i(x2 - xi)pi] exp[i(y2 - Z/Ops] exp[±M \ n - n \ ] (3.10) 
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where r = r'^Cr Ce and it is dimensionless. Further 

Ti = mm(ti,t2), ^3 = min{t3,t4). 

So far only 5-functions have been integrated. To integrate over remaining stochastic 
time variables, one needs to prove convergence of the integrals ( p.lO| ). This is easily done 
realizing that there exists a unitary transformation which diagonalizes the matrix (J) + m) 



( i\fp^ + m 

iy/p^ + m 











+ m 












\ 



U{p). (3.11) 



—i\fjP + 771/ 

After integrations over ti and t^, where it is necessary to integrate separately and then 



to add contributions from the regions > Ti and Ti > t^, one arrives at 



2r 



i2np (2vr)^ 



exp[z(a;2 - Xi)pi] exp[i(?/2 - yi)p3 



X 



X Tr|exp[-(2jf^^ + 2^^ + 4m) T] 



exp[(2;^^ + 2;^3 + 4m)T] 



exp[(2|^^ + 2m±M)T] 



{2^-^ + 2jf^3 + 4m) (2jf^3 + 2m T M) {2^^ + 2m ^ M) {2^-^ + 2m ± M) 
exp[(2j^^ + 2jf^3 + 4m)T] -1 exp[(2^ + 2m ± M)T] - 1 \] 



+ 



{2^^ + 2jf^3 + 4m) {2^^ + 2m ^ M) (2jf^3 + 2m ± M) {2^-^ + 2m ^ M) / J 

(3.12) 

where it has been put T = Ti = T3 in accordance with standard stochastic quantization 
procedure [0 . To carry out the limit T ^ 00 one has again to diagonalize the exponential 
operators by the unitary transformation ( p.ll|) . The result is 



{7P,r^f,ilj^) = 2r 



d^pi f d'^p3 



Tr 



{27rp (2vr) 
1 



4 exp[2(a;2 - xi)pi] exp[z(|/2 - 2/1)^3] 

1 1 M 



+ 



(^;^^+;^3 + 2m \^2|^3 + 2m=FM ' 2j(^^ + 2m =f M j J ' ^^'^^"^ 

where the signs =fM correspond to e^*^^ in ( p.22|) . In the case of correlations with e"'"^'^ 
it is necessary to assume 

2m > M (3.14) 

to assure convergence of ( p.l2| ) when T +00. 

The result ( p.l3|) is a four-point function and it is desirable to reduce it to a two-point 
function. To achieve this, new momenta have to be introduced 



p* = p^+p^ ^ p = p^-p.^ 
which substituted to ( p.l3| ) transform it to 



= 2r 



dY f d^p 



{2np (27r)^ 

exp [ip*[{x2 - xi) + (1/2 - yi)]] exp[iip[(x2 - Xi) - (^2 - Vi)]] 

1(1 1 M 



(3.15) 
(3.16) 



Tr 



+ 



p* + 2m \p* — p + 2m =p M p* + p + 2m =p M 



The space coordinate dependence splits into two parts. In one exponent one has the sum 
of relative coordinates of the fields ip and - this corresponds to p*. In the other exponent 
the "relative coordinate of relative coordinates" corresponds to p. It is useful to introduce 
also a new set of Euclidean space variables 

z = [{x2 - xi) - {y2 - yi)] 

u = [(x2 - xi) + (?/2 - yi)] (3.17) 
V = [yi + y2] 

W = [Xi + X2] 

The expression ( |3.16|) depends on z and u only. Variables v and w are chosen to complete 
the set of linearly independent variables. Alternative choices are 

v' = [x2 + xi + y2 + yi] , w' = w 

or 

v" = v' , w" = [y, + y2]. 

One possible way, how to reduce the number of space variables, is to integrate (|3.16|) 
over u. This yields 



{{ij^ri^pr)) ^ jdu i^^r^^r) = ^rj dYj ^ kpi e 

Tr{^- ( + ] ] . (3.18) 

and subsequent integration over p* leads to 

W„,rrp,,p')) = 2rf ^ (?^) . (3.19) 

J {2tt) V m / p^ + {2m ^ M)^ 

The further reductions of space variables imposed by conditions v = w = lead to 



where all coefficients were absorbed into the constant A^. By means of the same procedure 
applied to the conjugated field combinations one gets 

mx)r{-x)My)^'i-y))) = ^/ ''''''''' f + i2!nTMr ' ^'-''^ 

which is in accordance with the condition that the overall center of gravity is at rest. 



4 Long time limit 

In the previous section an expectation value and its limit when the stochastic time goes 
to infinity was calculated. To be able to ascribe some meaning to this result, one has to 
know that it is based on proper long time limit behavior of the whole process. This means 
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- one has to prove that the distribution of the solutions to the equations ( |3.1| ) has a hmit 
when the stochastic time goes to infinity and this hmit is unique. 

To achieve this goal the Fokker-Planck equation cannot be used, as it is derived under 
assumption that the relaxation time of fiuctuations is very small compared to observation 
time. The considered stochastic process of Grassmann numbers allows both signs in 
the auto-correlation function (e^*^"^) and thus it is desirable to keep full control of the 
time evolution rather then to use an approximation. A significant virtue of stochastic 
quantization is that it offers the whole arsenal of stochastic differential equations. It 
is possible to use the stochastic Liouville equation P, 0. The idea is to use the well 

known hydrodynamic equation P = — div p v in the space of solutions to (|3.1|) . The p 
becomes then the density of solutions ip and ip and there exists a theorem |^ saying 
that {pi^ipyi^yt)), i-e. the solution density averaged over all fluctuations r],^,f],^ is the 
probability distribution of ip, at stochastic time t. It is possible to prove the existence 
and uniqueness of the long time limit from the properties of such probability distribution. 
The Liouville equation in this case reads 

d d d 

— p{^, ij, t) = +^ [tp, p{ij, ij, t)] + — pitp, tp, t)], (4.1) 

where the plus sign on the right hand side comes from derivation of the Liouville equation 
with anticommuting variables 

Substitution of the right hand sides of ( ^.11 ) into ( |4.1| ) yields in p-representation 

dp \ . . , , (9 , , .yd 

+ {p + m)^f,i;^-— + {^ + 711)^^^)^- 



-Tn (e/. + V'.)^ - rn i^', + V'.)^^ P(^, V-, t) (4.2) 

where 77/ and f// are the Fourier pictures of stochastic terms and they have the same 
correlations as in (I^H) -(|2]2|). The probability distributions ( glSj^lOD have also to be 
transformed, e.g. the ( |2.18| ) transforms as follows 

V (e(p, t),ap', t), r]{p, t), nip', t)) = J^lVi^x, t),^{x', t), 7]{x, t), f]{x', t))] 
e^p-e-^'-' P{ax,t),ax,t),v{x,t),v{x,t))6{x-x') 



i2np (2vr) 
d'^x 



{2nrJ (2vr) 



P{ax,t),ax,t),v{x,t),n{x,t)) (4.3) 



e 



In order to solve (|4.2| ) it is necessary to transform it to the interaction representation. 
The stochastic time independent part of the operator on the right hand side of ( [4.2|) can 
be denoted as 

d - d 

Ao = -8m + + m),^^^-— + (^ + m),^V^,— (4.4) 

and the rapidly varying part as 

Olpu Olp^ 
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After transformation 

a{t) = e-*^V(t) , ^(0) = p(0) (4.6) 
the equation ( |4.2| ) turns into 

a (t) = e-*^°Ai(t)e*^« a{t) = -V{t) a{t) (4.7) 

where 

V{t) = e-*^« A,{t) e*^°. (4.8) 
Solution to ([4.7|) can be written as the time ordered exponential 

(T{t) = \exp{- f' V{t') dt'}]a (4.9) 
Jo 

which can be expanded as 

a{t) = a- dtiV{ti)a+ dtJ dt2V{ti)V{t2)a - dti dt^ dt^V {ti)V {t2)V {t^)a 
Jo Jo Jo Jo Jo Jo 

pi /'^S /'^S 

+ dti dt2 dh dUV{ti)V{t2)V{h)V{U)a- ... (4.10) 
JO Jo Jo Jo 

One gets rid of all the terms in ( |4.10| ) except of those of the order 4n(?T. = 0, 1, 2, ..) 
by calculation of the average {(T{t)) over Grassmann stochastic variables. The remaining 
terms are nonzero provided that they have been averaged in the same way as in (|3.7|) with 
the weight ( |3.8|) . The first two nonzero terms are 

{So) = {a) = {p{t = 0)) ^ 1 (4.11) 

which is in fact normalization of probability distribution and 

= f dti f dt2 f'dts [' dUV{ti)V{t2)V{t3)V{U)a. (4.12) 
Jo Jo Jo Jo 

256 terms of the product in 5*1 are reduced by integration over stochastic terms to the 
expression 



-'0 oi/jApi) dtl^Api) Jo 



_^tiAo^±M\ti-tz\ 



dilJy{pi) d^J^ipi) Jo dip^^^q^) dip^iqs 

+ 4r [dtie~^^^° - e*^^°/tit3e~*«^° - e^^^Oe^^^^'^^"*^' (4.13) 

Jo d^^{qi)d^f,{qi) Jo d^^{p3)dtpu{p3) 

where the time integrations have been carried out with the weight ( |3.8| ) and r = rj^Cr Ce 
is dimensionless as in the (p.lO|) . The absolute values of time difference in ( [4.13| ) come 
from the fact that although the time sequence is ordered, the Markovian process can run 
in both directions. To cope further with ( [4.13| ) it is necessary to commute the exponentials 
with the derivatives. This leads to 

{Si) = ir dti / dt3 —— exp[ti{:^^ + m ± (ixp[ti{:^-^ + m ± 

Jo Jo oipi, 1 I oip^ 

d M M d 

■ exp[t3(^3 + mT -Tr)]p^ exp[t3(^3 + T -rr)]^ 



d , M M,, (9 



r a IVl M 

+ Arjdt^j^dti — -exp[ti(/^^ + m =F — )]z.r exp[ti(|^^ + m =F y)]c^^ 



d M M d 

■ exp[t3(^3 + m± — )]p^ exp[t3(^3 + m ± — )], 



+ cc{p ^ q) (4.14) 
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where one has to integrate separately and add contributions from the two possible time 
orderings ti > and ts > ti. The cc part originates from the second term in ( [4.13|) . 
It is now necessary to calculate an integral of the type 

I'J, = fdr{e^na,{e^n.s. (4.15) 
^ 

By integration per partes one gets 

f^dT{e^^)^p{e^^),s = {e^la^.A~l{e^^),s\l - dT{e^-)^^ A'l {e^^),, (4.16) 
and provided that [A, -B] = 

dT{e^^)^,{e^^),, [5,p5,, + A-lB^s] = (e^-)„^ A^^^ (e^-),^ |* (4.17) 





which results in 



I f^r(e^O„^ (e^07<5 = + 5^^B,,]-' (e^^^^ (e^0.<5 lo (4-18) 

•J 



In order to calculate ( [4.14| ) one has to substitute A = B = /^^+m±^ (resp. 
A = B = ^3+m=i=^). At this moment it is useful put 1^-,^ = (jf^*+|^)/2 and ^3 = {jl)*—f)/2 



and to apply the procedure outlined in ( |3.15| - ^.191 ). This is equivalent to substitution 



= - ^3 = 'p/2. Application of ( ^181) to (|]T|) then leads to 



ft d i) M i> M d 

{Si) = 2r dti —— exp[ti{- + m ± —)]^r exp[ti{- + m ± —^^ 

Jo OWy II II 



exp[ti(-- + mT — )]p^exp[ti(-- + mT — )i 



+ "^^Idts^iltJ^u. + ?>u4^^ + 2{2mTM)b,r^^X^ 

i M i M \ d 

exp[t3(| + m T _)]^^exp[t3(| + m ^ —)\^^ + O^^^(exp(0))| ^ 

■ exp[t3(-77 + ± ^)]p7 exp[t3(-- + m ± — )], 



+ cc{p^q) (4.19) 

The terms (9(exp(0)) come from the lower limits of the integrals. It is only their behavior 
as functions of the stochastic time which is relevant for further calculations. Multiplication 
of matrices and subsequent time integration lead to 



{Si) =re 



di)y di)^ m (2m =f M) [p^ + (m =f ^)2] dipp dip^ 



+ Of:;(exp[t(2;^ + 2m±M)]) 



'''' d^u d^/j^ di/jp di/j^ 

d d d d 
dip^ dip^ dipp dij^ ' 



+ 0^^{exp[t{-2^ + 2m±M)])^^;j-^^;j^. (4.20) 
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This expression can be written in the form of the integral 

, ~„„. d d d d 
dtpu dip dipp dip„ 



=4m / drWi^^ir)^ — ^— (4.21) 



where 



WrAT)= re 



+ Oi^"p{[-2^ + 2m±M] exp[r(-2|^ + 2m ± M)]) 

+ 0^;([2|^ + 2m±M]exp[r(2;^ + 2m±M)]). (4.22) 

The equation 

(.)^l + 4™/;.rtT-M^^A^ (4.23) 

can be then regarded as truncated expansion of the solution to the equation 

d d d d 
' dip^ dip^ dipp dip„ 



(a) = AmWi^;{t)^~^~{a). (4.24) 



After transformation back from the interaction representation, the equation for (p) be- 
comes 

(P> = A„(p«)> + 4,n J_ e-. (4.25) 

and after commutation of the field derivatives with the operator e~^^° one obtains 

d d d d 
' dipy dip dipp dip^ 



ip) = Aoipit)) + 4mPy,7(t)^ — {pit)) (4.26) 



where 

^'-^W - - a^"(expM-2,„±M),). (4.27) 



To investigate the asymptotic properties of solutions to ([4.26| ) it is necessary to find 
its eigenstates and corresponding eigenvalues. To this purpose it is useful to introduce 

d 



(4.28) 

d 



s 



where U is the unitary transformation which diagonalizes the matrix and hence also W . 
The operators (|4.28| ) are anticommuting 

{a^at} = 5^ , {b.,btM = 5^. (4.29) 
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To convert ( f4.26| ) into algebraic equation one has to introduce further the states 



a*^ I 0) = , (0 I 



, 



0) = , (0 I b"^^ = 0. 



(4.30) 



The search for the eigenstates of the operator on the right hand side of ( [4.26|) leads to 
the equation 



= -8m + TW^a+a'^ + A^^b+^^b^ + 4m W^Ja'^b^a^'b^ 



R 



(4.31) 



where 
and 



To diagonalize ( [4.31|) it is necessary to introduce the similar transformation 



O = expf-W.^^^'b^b, 



This transformation changes the operator R on the right hand side of ( |4.31|) to 

O R = -8m + A<>^a^ + A^'b^'^b^. 
The kets are transformed as 



(4.32) 



(4.33) 



\P) = O \ P). 

The equation (|4.33| ) has 2^ different eigenstates with generally complex eigenvalues. 
The real parts of the complex eigenvalues are discreete equal to — m, —2m, ... — 7m, the 
imaginary parts are functions of the momentum p. The only two eigenstates with pure 
real eigenvalues and —8m are nondegenerate. 

The eigenstate corresponding to the eigenvalue Aq = can be written as 



^o) = n at n bt- I 0). 

u=0 /x=0 



(4.34) 



Thus any solution to ( [4.1[ ) averaged over stochastic terms in this representation (and this 
order of approximation) can be written as linear superposition of eigenstates 



255 



(p(t))) = j:e'^'c, \P^) 



(4.35) 



i=0 



and in the long time limit {t ^ oo) survives only the first term in the sum ( [4.35| ). 
Simultaneously vanishes the term (9^^(exp[t(— 2m ± M)]) in (^4.27 ). To untangle the 
transformations done, it is necessary to calculate 

3 3 

I Po) = I Po) = exp (W.^a'^b^a'^b,) fl aif[ h^^ \ 0). (4.36) 

a=0 /3=0 



The result is 



Po) = fip) exp [2 iW-Xl^^H^^'H] I 0). 



(4.37) 
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where 



Hp) 



[4p2 + (2m ^ Mf]^ + {2m T M)' 



(4.38) 



w? (2m T My \p^ + {m^ f )2]2 

The transformation back to nondiagonal ;^ matrices is straightforward. 

Thus in the long time hmit the expectation value of any functional {F[ip, t/S]) is given 

by 

(F[V^,V^]) = J F[^,V'] Po[ip,ip] Di^D,/; (4.39) 

or 

JF[iIj,'4j] exp{-r-^ Sl-^,^]) Dip Dilj 



(F[V',V']) 



where 



A 



S[^, ij] = ( -^,J^^ - i-p' + (2m T My]6.J,,)ij,ij^ . 



(4.40) 



(4.41) 



5 Conclusions 

Described method of the stochastic quantization of the fermion field yields as the expec- 
tation value of the combination -ip ^ip'^'ip pip^ the scalar boson propagator. The procedure 
leading to this result is based on the concept of Markovian processes of Grassmann vari- 
ables. The notion of derivative and integral in the Grassmann number calculus allows not 
only to introduce the concept of probability but also to extend it to stochastic processes. 
The exphcit form of one particular process has been constructed. This process is analog- 
ical to classical Ohrenstein-Uhlcnbeck process in the sense that all higher moments are 
products of the second moments. However, contrary to the classical process, in the case of 
Grassmann variables both negative and positive exponents in autocorrelation exponential 
are admissible. This sign appears in the resulting formulae as the sign of the term which 
can be interpreted as the mass defect. The negative mass defect stems from the process 
which has the "forbidden" sign. 



Calculation of the field correlations proceeds according to standard stochastic quan- 
tization procedure. The new problem arises with the necessity to reduce the number of 
dimensions both of the stochastic time space and of the Euclidean space. This is con- 
nected with transition from a four point function to a two point function. The reduction 
of two degrees of freedom in stochastic time follows quite naturally from the properties of 
the probability distribution of the stochastic terms. The manifold where it is nonzero is 
twodimensional. Staying on this manifold means, that the product ipa'4'°' is always taken 
in the same stochastic time instant, which is natural if one wants to interpret it as one 
particle composed of ip and ip. 

The initial system consisted of four fermions and the final system represents a boson 
field. Thus the coordinates of individual constituents disappear from description of the 
final system. However, the stochastic quantization provides no mechanism to do this. 
Hence the redundant degrees of freedom in Euclidean space are reduced by an ad hoc 
procedure which is based on transformation of the momenta to the center of mass system of 
'4'aip°' and its conjugate (which is identical in this case). The whole expression (ip^ip^ipf^ip^) 
is then integrated over space variable corresponding to the momentum of the overall center 
of mass p*. Consequently p* appears as an argument of a 5-function and it is set to zero 
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by integration. The whole system is not considered to be on mass shell. Two other space 
variables v and w ( p.l8|) are also set to zero which has the meaning of putting the pair %p 



and into one space point x resp. y. As a result the final boson propagator depends only 
on the distance of two points in Euclidean space and the systems (..) and (..) propagate 
with opposite momenta. The energy-momentum conservation constraint enters the whole 
description only in this sense. 

The dynamics of the fermion field affected by stochastic process ( p^.l8| , |2.igD is described 
in the long time limit by (|4.41| ) . It consists of the part which is up to the coefficient of the 
mass term equal to free scalar boson Lagrangian and the term which is analog to the scalar 
field term in Stueckelberg's Lagrangian of a massive vector boson field. The functional 
( [4.39| ) does not act on any arbitrary scalar function but only on spinor functions. Hence 
this term will lead in case of a scalar constructed from spinor functions {F[^ ^ip'^if) ^ip^^]) 
to Tr{^) and ( [4. 411) will turn to free scalar boson Lagrangian with no selfinteraction. In 



case of other functions constructed from spinors, ( [4.40[ ) gives generally nonzero results - 
e.g. Fl'i/j'j'^ilj tp'j'^ilj] results in constant —5^'^/m{2m =p M). However, it has to be stressed 
that the probability distribution of stochastic terms is just the most simple one satisfying 
conditions of a Markovian process with no ambition to describe any particular dynamics 
in the long time limit. 

The long time behaviour of the whole system is studied starting from the expansion of 
the solution to the Liuville equation. Thus the quality of approximation by Fokker-Planck 
equation can be controlled. The convergence of the expansion is given by numerical values 
of r and fermion mass m as the coefficients corresponding to m (2m=l=M) in denominator of 
5*1 (|]2g) become equal to (6m=FM)(2m=FM), (10m=FM)(6m=FM)(2m=FM), ..etc 
in higher terms 5*2, S^, ... The Gaussian behaviour of stochastic terms implies that higher 
order terms should reduce to products of free "boson Lagrangians" (each containing the 
term discussed above) with masses increasing by 4m. However, then the procedure leading 
to reduction of the degrees of freedom has to be generalized and the result will depend 
on it. 



The whole calculation closely follows the standard stochastic quantization procedure. 
The only two differences in the approach consist in different statistical properties of the 
stochastic terms and in reduction of the degrees of freedom in final formulae. The later 
difference is clearly connected with desired "coarse graining" of the description. On 
the other hand the statistical properties of stochastic terms in Langevin equations are 
responsible for the fact that in the long (stochastic) time limit the free fermions disappear 
and the combination {ifj^ip'^'il) f^xjj^) propagates as a scalar boson field. Thus the stochastic 
terms can be regarded as the external field coupling the fermions together. The question 
whether this interaction is an effective result of a standard interaction via gauge fields lays 
beyond the scope of the described procedure as the statistical nature of stochastic terms is 
imposed ad hoc. Any hint that the given method refiects real dynamics would generate a 
puzzling problem - mere statistical properties of stochastic term make difference between 
standard stochastic quantization of free fermion field and description of interacting system 
which gets quantized " along the way" in course of calculation. 

It is a pleasure to thank Jiri Rames and Jiri Chyla for reading the manuscript and 
also to acknowledge very stimulating discussion with Jiri Hofejsf. 
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